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ABSTRACT 

Input  Control  plays  a  critical  role  in  regulating  the  release  of  jobs  into  a  production  system. 
For  a  just-in-time  system,  we  model  the  input  control  decision  as  a  bicriterion  problem 
which  requires  the  maximization  of  the  sum  of  job  release  times  subject  to  minimum  total 
job  tardiness. 

The  proposed  solution  method  uses  a  constructive  approach  in  which  jobs  are  grouped 
into  blocks  for  simultaneous  processing.  The  membership  of  each  block  is  updated  in  an  iter- 
ative manner  by  solving  a  sequence  of  subproblems.  Our  computational  experience  indicates 
that  this  approach  yields  substantial  improvement  in  job  release  times  without  affecting  the 
tardiness  values  adversely.  This  result  is  of  interest  to  an  operating  manager  who  is  responsi- 
ble for  controlling  work-in-process  levels  and  earliness  costs  in  addition  to  meeting  due  dates 
effectively. 


1      Introduction 

In  many  manufacturing  systems,  the  release  of  jobs  into  the  production  shop  floor  is  regulated 
by  means  of  an  input  control  mechanism.  The  purpose  of  this  mechanism  is  to  allow  jobs 
to  be  introduced  selectively  into  the  system  instead  of  permitting  all  of  them  in  as  and 
when  they  become  available.  In  most  materials  requirements  planning  (MRP)  systems,  this 
mechanism  results  in  important  savings.  First,  the  work-in-process  levels  are  reduced  leading 
to  less  congestion  and  confusion  on  the  shop  floor.  Second,  by  deferring  jobs  required  at  a 
future  date,  the  system  is  less  susceptible  to  demand  variations,  and  is  therefore,  less  nervous. 
Third,  it  simplifies  scheduling  and  smooths  the  output  obtained  from  the  system. 

Input  control  is  also  imperative  in  many  flexible  manufacturing  systems  (FMS)  because 
of  the  limited  buffer  space  available.  Input  control  may  also  be  an  integral  part  of  the 
hierarchical  planning  and  scheduling  of  an  FMS.  Other  instances  when  input  control  is 
required  are  for  shop  load  leveling  (Irastorza  and  Deanne  1974),  load  balancing  (Shimoyashiro 
et  al.  1984),  part  and  tool  loading  in  an  FMS  (Stecke  1983),  etc. 

While  there  are  may  practical  benefits  of  using  input  control,  Baker  (1984a)  notes  some  of 
the  possible  drawbacks  as  well.  By  removing  some  of  the  available  options  to  the  scheduler, 
input  control  may  be  counterproductive.  Specifically,  from  a  scheduling  perspective,  if  the 
scheduling  objectives  relate  to  regular  measures,  i.e.,  if  the  schedule  costs  are  nondecreasing 
in  job  completion  times,  then  delaying  job  input  could  easily  result  in  suboptimal  solutions. 

The  just-in-time  (JIT)  approach,  however,  promotes  the  notion  that  jobs  should  be  re- 
leased as  close  to  their  due  dates  as  possible  even  if  they  are  available  for  scheduling  much 
earlier.  Many  scheduling  objectives  relevant  to  JIT  systems,  such  as  minimizing  earliness, 
are  nonregular.  Given  a  set  of  jobs  available  for  scheduling  at  any  point  in  time,  the  basic 
approach  adopted  in  such  systems  is  to  delay  the  release  of  noncritical  jobs  as  much  as  possi- 
ble while  insuring  at  the  same  time  that  critical  jobs  are  taken  up  without  delay.  Therefore, 
given  a  set  of  jobs  available  for  scheduling,  the  problem  facing  the  production  planner  is  to 
identify  the  critical  jobs  which  need  to  be  input  immediately.  In  addition,  (s)he  needs  to 
determine  the  future  release  times  for  the  deferred  jobs. 

In  this  paper,  we  model  this  decision  as  the  problem  of  maximizing  the  release  times 


of  nonurgent  jobs  subject  to  the  minimization  of  the  sum  of  tardiness  of  urgent  jobs.  This 
bicriterion  problem  considers  two  objectives  ordered  hierarchically.  The  primary  objective  is 
to  minimize  the  sum  of  tardiness  of  urgent  jobs.  Subject  to  this  requirement,  the  secondary 
objective  deals  with  the  maximization  of  the  release  times  of  nonurgent  jobs.  Note  that 
while  the  primary  objective  is  a  regular  measure  of  performance,  the  secondary  objective  is 
nonregular. 

While  this  is  the  first  known  attempt  to  combine  these  two  objectives,  the  literature 
on  the  job  shop  tardiness  problem  considered  independently  is  quite  rich.  In  view  of  the 
exponential  complexity  of  this  problem  (see,  for  example,  Rinnooy  Kan  1976),  prior  research 
on  this  problem  has  primarily  considered  the  use  of  priority  dispatching  rules.  [See  Baker 
1984b  and  Vepsalainen  and  Morton  1987  for  an  extensive  discussion  on  this  research.]  In  a 
recent  study,  Raman,  Talbot  and  Rachamadugu  (1989)  developed  an  implicit  enumeration 
approach  as  well  as  a  decomposition  based  heuristic  for  solving  this  problem. 

While  the  problem  of  maximizing  the  sum  of  job  release  times  in  a  job  shop  has  not  been 
addressed  directly,  there  is  considerable  literature  on  this  subject  for  single  machine  systems. 
Note  that  for  a  single  machine  system,  this  objective  reduces  to  that  of  minimizing  earliness. 
This  problem  is  shown  to  be  NP-complete  by  Chand  and  Schneeberger  (1986).  Baker  and 
Scudder  (1990)  provide  a  survey  of  the  work  done  on  this  objective  and  its  variants.  Ahmadi 
and  Bagchi  (1987)  have  extended  this  objective  to  the  case  of  a  two-machine  flow  shop.  This 
paper  is  organized  as  follows.  In  §2,  we  present  an  integer  programming  formulation  of  the 
input  control  problem  and  discuss  its  decomposition.  The  solution  procedure  is  discussed  in 
§3.  Our  computational  experience  with  this  approach  is  described  in  §4.  We  conclude  in  §5 
with  a  discussion  of  the  experimental  results.  The  notation  used  in  this  paper  is  given  in 
Appendix  1. 

2      The  Input  Control  Problem 

The  input  control  problem  can  be  formulated  as  below. 
ICP 


n 
Maximize  z2  =  ^fj  (1) 

3=1 

subject  to 

N 

y^Tj  =  zi, where  Z\  —  mina€s{T{cr)}  (2) 

3=1 

xUv  -  x3i  +  Af (1  -  <^tuv)  -  put,  >  0,   V;',  u€  J,   and  (j,  z),  (u,  v)  G  ^m,   m  =  1, . . . ,  M    (3) 
Zji  -  zUv  +  M6jiuv  -  pj{  >  0,   Vj,  u£j,   and  (j,  i),  (u,  u)  G  //m,   m  =  1, . . . ,  M       (4) 

f   I*;      if  /  =  1 
a?iz  -  Pj/  >  \  .  .  (5) 

Xj,    otherwise;  Vj,  i,  and(z,  /)  E  Cj 

XjiVj  +  Ej  -  T3  =  d3,Vj  (6) 

Sjiuv  G  {0, 1 },  Vj,  i,  u,  u;  Xji  >  OVj,  i;   r;,  £;,  7j  >  0,  Vj  (7) 

where  Xji  is  the  completion  time  of  operation  i  in  job  j,  and  <5JJUV  equals  1  if  operation  i  in 
job  j  precedes  operation  v  in  job  u,  and  is  zero,  otherwise. 

In  this  formulation,  Equation  (1)  refers  to  the  secondary  objective  of  maximizing  the 
sum  of  job  release  times.  Constraint  (2)  specify  the  primary  objective  of  minimizing  total 
tardiness.  In  the  given  expression,  T(a)  refers  to  the  tardiness  incurred  in  schedule  <r,  and  S 
refers  to  the  set  of  all  schedules.  Constraints  (3)  and  (4)  are  the  disjunctive  constraints  that 
ensure  that  no  more  than  one  operation  is  processed  on  a  machine  at  any  time.  Constraints 

(5)  ensure  that  any  operation  can  be  scheduled  only  after  its  predecessor  operation  is  com- 
pleted, and  the  first  operation  of  any  job  is  started  only  after  the  job  is  released.  Constraints 

(6)  define  tardiness  while  Constraints  (7)  specifies  the  nature  of  the  variables. 

ICP  generalizes  two  well-known  problems.  The  primary  objective  addresses  the  mini- 
mization of  mean  tardiness  in  a  job  shop.  If  it  is  feasible  to  complete  all  jobs  by  their  due 
dates,  then  the  problem  reduces  to  maximizing  the  sum  of  job  release  times  subject  to  all 
jobs  being  completed  by  their  due  dates.  Maximizing  Y,)€jrj  xs  equivalent  to  minimizing 
EjeJ-K  -rj).  Also, 

(dj  -rj)  =  J2(XH  ~  x3,i-i)  +  Pji- 

t=2 

If  we  denote  the  waiting  time  of  job  j  after  operation  i  by  wJt,  then 

xJt  —  av,,_i  =  v)j,i-i  +  Pji- 


Hence, 

(dj  -  rj)  =  {d,  -  Cj)  +  ]T  Wji  +  Y,pji 

t=i  «=i 

The  last  term  on  the  right  hand  side  is  a  constant.  Hence,  minimizing  J2j€j(dj  -  rj)  is 
equivalent  to  minimizing  T,j€j(dj  -  Cj)  +  E^ESi"1  u>i»-  Tnis  objective  has  been  studied 
by  Ahmadi  and  Bagchi  (1987)  for  a  two  machine  flow  shop.  In  the  case  of  a  single  machine 
job  shop,  this  objective  reduces  to  minimizing  total  earliness  which  is  a  well  researched 
problem  (see,  for  example,  Baker  and  Scudder  1990). 

Both  job  shop  mean  tardiness  problem  (see,  for  example,  Rinnooy  Kan  1976)  and  the  sin- 
gle machine  earliness  problem  (Chand  and  Schneeberger  1986)  have  been  individually  shown 
to  be  NP-complete.  Consequently,  we  need  to  use  heuristic  methods  for  solving  any  problem 
of  reasonable  size.  The  proposed  approach  is  based  on  identifying  a  sequence  of  subproblems 
within  ICP,  and  constructing  a  schedule  by  solving  these  subproblems  iteratively. 

3      Embedded  Subproblems 

In  the  rest  of  the  paper,  we  will  assume  that  all  jobs  are  numbered  in  the  order  of  increasing 
due  dates,  i.  e.,  if  v  >  u,  then  dv  >  du.  Intuitively,  it  is  clear  that  jobs  that  are  either 
already  tardy  or  are  in  the  imminent  danger  of  being  so  need  to  be  released  into  the  system 
immediately,  while  the  input  of  nonurgent  jobs  can  be  deferred.  Consider  such  a  partition 
of  J  into  the  set  of  urgent  jobs  J\  and  the  set  of  nonurgent  jobs  J2,  J1UJ2  =  J,  an^ 
Jx  f]J2  =  0;  and  a  schedule  a  =  {JUJ2)  in  which  all  jobs  in  Jx  are  completed  before  any 
job  in  Ji  is  started,  i.  e., 

max]€Jx{c3}  <  minj€j2{rj] 

Let  c(.)  denote  the  completion  time  of  jobs  in  a  given  set  in  schedule  a.  Then, 
Remark  1   In  an  optimal  solution,  ifJ2jejTj  >  0,  then 

11)  Tj  <  c{Jx\j)  +  Pj  -  dj,Vj  e  Jx,   and 
Hi)  J\  starts  at  time  zero. 


PROOF:  Refer  to  the  Appendix  2. 

For  a  given  partition  (J\,J-i)  that  satisfies  the  conditions  given  in  Remark  1,  ICP  sep- 
arates into  two  problems,  ICP1  that  considers  only  J\  and  ICP2  that  considers  only  J2  as 
given  below. 


ICPl(Ji) 


subject  to 
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Maximize  z2  =  2_.  rj 
jeJi 


(3)  -  (7),  Vj  G  Jr. 


(8) 


(9) 


ICP2(J2) 


subject  to 


and, 


Maximize  z2  =  tJ  r? 


T,  =  0,  Vj  G  J2 
7\,  -  maij€j,{cj}  >  0  Vj  G  J2 

(3)  -  (7),  V;  G  J3. 


(10) 

(11) 
(12) 


Constraints  (11)  and  (12)  insure  that  the  partition  (J\,J2)  is  feasible.  Note  that  the 
dual  objectives  apply  only  to  ICP1,  while  ICP2  considers  only  the  objective  of  maximizing 
the  sum  of  job  release  times. 

In  general,  any  schedule  a  optimal  to  ICP  will  have  inserted  idle  times  because  the 
secondary  objective  is  a  non-regular  measure  of  performance.  This  results  in  job  blocks  as 
shown  in  Figure  1 .  Two  jobs  ji  and  jn  belong  to  the  same  block  if  a  consists  of  a  subsequence 
Ui,J2,-  •  -,jn)  such  that 

cn  >  Gj,_,  >  rjn   I  =  2.  ...,n 


It  follows  that  for  two  adjacent  blocks,  B[k]  and  B[fc+i], 

maxi€Blk]{cj}  <  minj€B[k+11{rj}. 

Note  that  if  J\  is  not  empty,  then  it  constitutes  a  block. 

ICP2  can  now  be  reformulated  in  terms  of  blocks  as  follows: 


Maximize  z2{ICP2)  =  £    £   r3  (13) 

k=ljeB[k] 


subject  to 


J7 

^Vjk  =  1,   Vj  e  J2 
k=i 


maxjeB[k]{Cj}  -  rnin^B^irj)  <  0,   VA; 
m?"j€S[1]{rJ}  -  ma-r^jc^}, 
JUv  -  Xji  +  M(l  -  6jiuv)  -  puv  >  0,  Vj,u£Bk,  and  (j,i),(u,  v)  €  //m,  m  =  l,...,M 
Xjj  -  ;ruv  +  MSjiuv  -  pJt  >  0,   Vj,  u  £  Bk,   and  (j, »),  (u,  u)  6  //m,   m  =  1, . . . ,  M 

r3      if7  =  1 

Xj,    otherwise;  Vj,  i,  and(z, /)  G  Oj 


14) 

15) 
16) 
17) 
18) 

19) 


xjN]  +  Ej  -  T,  =  dj ,  Vj  e  Bk,  Vfc  (20) 

*i««t»  £  {0, 1},  Vj,z',w,v;  Xj,  >  0,  Vj,i;  rj,Ej,Tj  >  0,  V;  (21) 

For  a  given  partition  of  J2  into  A'  blocks,  {Bk}k-n  that  satisfies  (14)— (16),  ICP2  decom- 
poses into  K  independent  subproblems;  each  subproblem  addresses,  independently  for  each 
block,  the  maximization  of  Ylj€Bk  rj  subject  to  (17)— (21 ).  We  shall,  henceforth,  refer  to  this 
problem  as  the  block  maximum  release  time  problem  (BMRP).  For  a  given  sequence  o~k  of 
jobs  in  block  Bk-,  we  define  the  block  release  time  Rk,  block  completion  time  Cjt,  block  due 
date  Dk  and  block  processing  time  Pk  as: 

Rk  =  minjeBk{rj} 

Ck  =  maxj€Bk{cj} 

Dk  =  rnirij€Bk{dj\cj  =  Ck] 

Pk  =  maxjeBk{c:j}  -  minjeBk{sj} 


Clearly,  for  the  partition  and  the  schedule  of  jobs  in  each  block  to  be  feasible,  constraints 
(14)— (16)  need  to  be  satisfied,  where  (15)  can  now  be  restated  as: 

C[k]  -  %+1]  <  0,   VA; 

In  addition,  the  following  results  can  be  established  for  an  optimal  schedule. 

Remark  2  In  each  block  Bk,  k  =  1, . . .  ,K  in  an  optimal  schedule,  there  exists  one  critical 
job  j  such  that 

Cj  =  dj. 

PROOF:  Refer  to  the  Appendix  2. 

We  will  denote  the  completion  time  and  due  date  of  the  critical  job  in  block  Bk  by  Ck» 
and  Dk*,  respectively. 

4      Solution  Approach 

A  heuristic  solution  procedure  can  be  constructed  utilizing  the  results  obtained  in  the  previ- 
ous section  in  the  following  manner.  First,  we  assume  that  all  jobs  can  be  finished  on  time, 
and  set  J2  =  J\J\  =  0.  Next,  a  set  of  initial  blocks  is  constructed  through  a  maximal 
partition  of  J2,  and  the  maximum  release  time  problem  is  solved  for  each  block.  If  there  is 
no  feasible  solution  for  a  given  block  Bk  (which  implies  that  at  least  one  job  in  Bk  cannot 
be  completed  by  its  due  date  in  this  schedule),  all  jobs  in  Bk  are  transferred  to  set  J\,  and 
J2  is  updated  accordingly.  At  the  end  of  this  step,  a  tentative  partition  of  J  into  J\  and 
J~2,  and  a  partition  of  J2  into  blocks  that,  considered  independently,  consist  of  early  jobs  is 
obtained. 

Next,  problem  ICP1  is  solved  for  jobs  in  J\.  Subsequently,  we  solve  a  relaxation  of  ICP2 
obtained  by  ignoring  (15)  and  (16)  with  respect  to  the  blocks  in  J2  by  considering  them  in 
the  increasing  order  of  their  indexes,  and  sequencing  them  such  that  Ck»  =  Dk»  for  each  k. 
If  this  sequence  satisfies  (15)  and  (16),  then  from  Remark  4,  it  is  optimal  to  ICP2,  and  the 
algorithm  terminates.  Otherwise,  if  (15)  is  violated  by  (say)  block  Bk+i.k  >  1,  then  it  is 
merged  with  block  Bk-    This  requires  solving  BMRP  for  Bk  with  the  enlarged  set  of  jobs, 


and  updating  the  values  of  Ck*  and  Dk> .  Block  Bk  is  sequenced  such  that  CV  =  Dk* .  This 
exercise  is  repeated  for  all  such  blocks  that  violate  (15). 

On  the  other  hand,  if  (16)  is  violated,  then  jobs  in  B\  are  transferred  to  J\,  and  ICP1  is 
re-solved  with  the  respect  to  the  new  set  of  jobs  in  J\.  The  indexes  of  blocks  in  Ji  is  updated 
accordingly,  and  the  process  is  repeated  until  (16)  is  satisfied  with  the  revised  makespan  of 

This  procedure  will  terminate  within  \J\  —  1  steps  yielding  a  sequence  which  provides 
a  feasible  partition  of  J  into  J\  and  J7-2,  and  also  insures  that  the  blocks  in  Ji  satisfy  the 
conditions  stated  in  Remark  4.  A  formal  statement  of  the  algorithm  is  given  in  Appendix  3. 
The  major  steps  in  the  algorithm  are  now  discussed. 

4.1  Determination  of  Initial  Blocks 

The  initial  set  of  blocks  is  obtained  from  a  maximal  partition  of  J .  A  partition  is  called 
maximal  if 

1.  each  block  in  the  partition  satisfies  the  following  condition:  If  u  and  v  are,  respectively, 
the  smallest  and  the  largest  job  indexes  in  the  block,  then  dv  —  pv  <  du,  and 

2.  minjeBk {d3  —  Pj]  >  maXj^B,,^  {dj},   k  =  2, . . . ,  K 

It  can  be  shown  that  in  an  optimal  solution,  if  the  job  that  is  completed  the  last  in  each 
block  has  the  largest  due  date  among  all  jobs  in  that  block,  then  the  resulting  partition  is 
maximal.  An  0(N2)  algorithm  for  obtaining  such  a  partition  is  given  in  Appendix  3. 

4.2  The  Block  Maximum  Release  Time  Problem 

For  each  block,  we  next  solve  the  block  maximum  release  time  problem  (BMRP)  -  finding 
the  sequence  of  jobs  in  each  block  that  maximizes  the  sum  of  job  release  times  subject  to 
the  requirement  that  they  are  completed  by  their  due  dates. 


4.2.1      Initial  Solution 

The  solution  method  proposed  for  BMRP  is  an  improvement  heuristic  for  which  an  initial 
solution  is  obtained  by  constructing  a  schedule  for  an  equivalent  problem  -  the  completion 
time  problem  (CTP).  The  objective  of  CTP  is  to  minimize  the  sum  of  completion  times  of 
jobs  subject  to  ready  times;  the  ready  time  of  job  j  is  given  by  r'i  =  dmax  —  dj,  where  dmax 
is  the  maximum  of  the  due  dates  of  jobs  within  the  given  block. 

To  see  the  equivalence  between  these  two  problems  consider  the  problem  shown  in  Figure 
1.  In  this  3-job  example,  we  have  d3  <  d2  <  d\  =  dmax  for  the  BMRP.  The  corresponding 
instance  of  CTP  is  constructed  by  defining  the  following  job  ready  times: 

Tj  =  0,  r2  =  «i  — -  «2j   and    r3  =  d\  —  ds 


INSERT  FIGURE  1  ABOUT  HERE 


Denoting  the  variables  in  CTP  with  primes,  it  can  be  seen  that, 

^max  —  Cmaxi     and    Cj   —  O-max        Tj, 

Therefore,  maximizing  YL3  rj  in  BMRP  is  equivalent  to  minimizing  J2j  c'j  m  CTP.  BMRP  can 
then  be  solved  by  first  solving  CTP  to  determine  job  completion  times  c^c^,  and  c'3,  and 
then  using  the  above  relationship  to  obtain  the  ready  times  r3. 

While  CTP  remains  a  hard  problem  even  for  a  single  machine  system  (Hariri  and  Potts 
1983),  a  feasible  solution  to  it  can  be  constructed  quite  easily.  We  generate  a  non-delay 
schedule  by  considering  only  those  jobs  which  are  available  at  a  given  machine  when  the 
scheduling  decision  is  to  be  made  at  that  machine.  Ties  between  competing  jobs  are  broken 
in  favor  of  the  job  with  the  shortest  imminent  operation  time.  The  order  of  operations 
generated  for  CTP  is  reversed  to  obtain  the  desired  sequence  of  jobs  within  the  batch  for 
BMRP.  This  sequence  determines  the  start  times  sJt  of  each  operation  i  within  each  job  j 
for  the  given  block  at  the  appropriate  machine. 

4.2.2      Schedule  Improvement 

During  the  improvement  phase,  we  consider  one  machine  at  a  time,  in  order,  starting  with 
the  bottleneck  machine.  For  each  machine,  an  attempt  is  made  to  reschedule  operations,  and 


the  revised  schedule  is  then  used  to  update  the  parameters  considered  for  the  next  machine. 
Thus,  this  phase  is  iterative;  it  terminates  when  no  further  improvement  is  possible.  The 
problem  to  be  solved  at  each  machine  involves  maximizing  the  sum  of  release  times  of 
operations  subject  to  their  completion  by  the  specified  due  dates.  Because  this  is  a  single- 
machine  problem,  maximizing  the  sum  of  release  times  is  equivalent  to  minimizing  total 
earliness.  [We  will,  henceforth,  refer  to  this  problem  as  the  earliness  problem  subject  to 
release  times  (ERP).] 

Note  that  an  operation  i  in  job  j  can  not  be  taken  up  before  its  ready  time  rJt  which  is 

given  by, 

{cJ)t_i     if  i  >  1 
0  otherwise 

Furthermore,  because  we  consider  only  one  machine  at  a  time  and  keep  the  sequences  at 
other  machines  unchanged,  the  due  date  of  operation  dji  is  determined  by 


d»  =  < 


Sjti+i     if  i<  Nj 


dj  otherwise 

Let  /  be  the  total  number  of  operations  on  any  given  machine.  For  the  ease  of  presenta- 
tion, we  introduce  the  following  notation.  Let  a,,  bx,  £,-,  /,  and  e;  be,  respectively,  the  ready 
time,  due  date,  processing  time,  completion  time  and  earliness  of  operation  i.  ERP  can  now 

be  stated  as: 

/ 

Minimize  v  =  y^  e, 

subject  to, 

fi-U>  at,   Vi 

et  =  6,  -  fi 

fi  -  h  +  M(l  -  6hi)  -  U  >  0,  Vz,  h,i£h 

h  -  fi  +  A/(l  -  Shl)  -  U  >  0,   V*\M  ±  h 

fi,a  >  o,  Vi 

where,  Shi  equals  1  if  operation  h  precedes  i,  and  is  zero,  otherwise. 
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Given  the  complexity  of  the  regular  earliness  problem,  it  is  unlikely  that  a  polynomial  time 
exact  solution  algorithm  is  available  for  ERP.  We  construct  a  heuristic  solution  algorithm 
which  is  a  modification  of  the  approximation  method  suggested  by  Chand  and  Schneeberger 
(1985)  for  the  regular  earliness  problem.  For  the  sake  of  completeness,  we  give  an  outline  of 
the  Chand  and  Schneeberger  heuristic  (CH)  below. 

The  algorithm  builds  the  schedule  backwards  by  using  the  Smith  heuristic,  i.  e.,  the  iter- 
ative step  involves  selecting,  from  among  all  operations  feasible  for  assignment,  the  operation 
with  the  smallest  tt.  The  schedule  developed  in  this  manner  is  optimal  if 

t[i\  >  fp+i],   i  =  1,...,/—  1. 

Suppose  cr1  =  {h&}  and  a2  =  {/i,i,cr}  are  two  subsequences  that  satisfy  the  above 
optimality  condition.  However,  at  the  next  step  the  operation  selected  is  g  that  yields  the 
subsequence  a3  =  {g,h,i,cr}  that  violates  this  condition.  The  heuristic  then  backtracks  to 
the  partial  solution  a1  to  consider  two  alternative  completions.  In  the  first  sequence,  g  is 
assigned  to  the  position  just  ahead  of  i  and  the  rest  of  the  schedule  is  developed  by  using  the 
Smith  method.  In  the  second  completion,  g  is  assigned  to  the  next  to  last  position  among 
the  unscheduled  jobs  and,  subject  to  this  assignment,  the  best  sequence  obtained  by  using 
the  Smith  method  is  found.  The  better  of  these  two  completions  is  selected. 

When  all  operations  are  available  at  time  zero,  checking  if  any  operation  is  feasible  for 
assignment  at  a  given  position  requires  only  insuring  that  /,■  <  6,.  However,  with  non- 
zero ready  times,  the  algorithm  needs  to  additionally  insure  that  an  operation  does  not  start 
before  its  ready  time.  More  importantly  from  the  computational  perspective,  the  assignment 
of  an  operation  to  a  given  position  should  insure  that  a  feasible  completion  is  available  for 
the  remaining  subproblem,  i.  e.,  there  is  at  least  one  schedule  in  which 

fi  -U>  at 

for  each  unscheduled  operation  i. 

Determining  whether  such  a  completion  exists  for  a  given  partial  sequence  is  equivalent 
to  solving  the  problem  of  minimizing  the  maximum  lateness  of  the  unscheduled  operations 
subject  to  ready  times  (I/l/rt  >  0/ Lmax).  In  particular,  a  feasible  completion  exists  if  and 
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only  if  the  optimal  Lmax  <  0.  Although  7/1/r,-  >  0/ Lmax  is  NP-hard  in  the  strong  sense, 
effective  algorithms  for  solving  it  have  been  proposed  by  McMahon  and  Florian  (1975)  and 
Carlier  (1982)  that  are  based  on  implicit  enumeration.  We  modify  Carlier's  algorithm  in 
order  to  improve  its  efficiency  in  the  context  of  ERP.  This  is  done  by  identifying  conditions 
that  lead  to  infeasible  completions.  Clearly  a  positive  lower  bound  to  the  overall  problem 
gives  one  such  condition.  Our  computational  experience  revealed  that  preprocessing  to  insure 
that  the  following  two  criteria  are  satisfied  also  improved  the  overall  solution  time: 

1.  t  >  a,  +  £,-,   for  all  unscheduled  i,  and 

where  t  is  the  start  time  of  the  most  recently  scheduled  operation. 

Modifying  the  start  times  and  the  completion  times  of  an  operation  of  any  job  on  a 
given  machine  results  in  revising  the  ready  time  for  the  succeeding  operation  and  the  due 
date  for  the  preceding  operation  for  the  same  job.  Thus,  in  general,  rescheduling  a  machine 
leads  to  updating  the  parameters  for  the  next  machine  to  be  considered.  Starting  with  the 
bottleneck  machine,  we  reschedule  each  machine  in  the  order  of  their  workloads  in  each  cycle. 
The  cycles  are  repeated  until  no  further  improvement  is  observed.  Although  it  is  difficult  to 
theoretically  guarantee  convergence  of  the  objective  function  value,  we  found  it  to  be  so  in 
our  computational  experience;  indeed,  the  bulk  of  the  improvement  was  obtained  within  the 
first  three  cycles. 

4.2.3      The  ICP1  Problem 

As  discussed  earlier,  ICP1  is  bicriterion  problem  with  the  primary  objective  of  minimizing 
mean  tardiness  of  all  jobs  in  J\.  Our  solution  method  consists  of  solving  the  mean  tardiness 
problem  first.  The  job  due  dates  are  then  revised,  if  necessary,  based  on  the  resulting 
schedule.  Given  these  due  dates,  the  maximum  release  time  problem  is  next  solved  following 
the  procedure  discussed  in  the  previous  section.  Note  that  such  a  decomposition  of  the  two  I 
objectives  can  result  in  a  suboptimal  solution  if  there  are  alternative  optimal  solutions  to 
the  mean  tardiness  problem. 
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In  view  of  the  strong  NP-completeness  of  the  mean  tardiness  problem,  we  propose  using 
the  heuristic  method  developed  by  Raman,  Talbot  and  Rachamadugu  (1989)  which  is  briefly 
described  now.  This  solution  procedure  is  a  schedule  improvement  procedure  for  which  the 
initial  solution  is  generated  by  applying  the  Modified  Operation  Due  Date  (See,  for  example, 
Baker  1984b)  rule  with  operation  due  dates  (ODDs)  set  loosely  at  the  maximum  values  that 
they  can  assume  without  delaying  the  corresponding  jobs.  Each  machine  is  then  considered 
in  the  order  of  their  relative  workload,  and  an  attempt  is  made  to  revise  the  schedule  of 
operations  on  that  machine  by  modifying  their  operation  due  dates.  Jobs  processed  on  all 
machines  are  ranked  in  the  nonincreasing  order  of  their  tardiness.  For  any  operation  in  a 
given  job  with  positive  tardiness,  first  we  determine  the  appropriate  interval  for  searching  for 
the  ODD.  For  each  possible  ODD  value  in  this  interval,  the  entire  system  is  rescheduled.  The 
value  which  yields  the  minimum  total  tardiness  is  returned  as  the  ODD  for  that  operation. 
This  step  is  repeated  for  all  other  operations  of  that  job  processed  on  the  machine  under 
consideration,  for  all  other  tardy  jobs  on  that  machine  following  their  rank  order,  and  for  all 
machines  in  the  system. 

Given  the  minimum  mean  tardiness  schedule,  the  revised  job  due  dates  are  given  by 

dj  <—  max(dj,Cj). 

Subject  to  these  due  dates,  we  solve  BMRP  for  the  jobs  in  J\.  Because  we  already  have  a 
feasible  schedule  initially,  we  can  go  directly  to  the  schedule  improvement  phase  discussed 
in  Section  4.2.2. 

Note  that  solving  the  mean  tardiness  problem  can  result  in  zero  tardiness.  In  such 
instances,  it  may  be  possible  to  defer  the  release  time  of  jobs  in  J\. 

5      Computational  Experience 

Two  sets  of  experiments  were  conducted.  The  first  set  evaluated  the  improvement  obtained 
in  the  sum  of  job  release  times  by  solving  the  dual-objective  input  control  problem  over  the 
single-objective  mean  tardiness  problem.  While  it  is  clear  that  explicit  consideration  of  this 
objective  should  result  in  better  solutions,  the  purpose  of  this  set  of  experiments  was  to 
examine  the  margin  of  possible  improvement  achieved  by  using  the  suggested  algorithm. 
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The  objective  of  the  second  set  of  experiments  was  to  judge  the  effectiveness  of  the 
proposed  solution  method  with  respect  to  known  upper  bounds  for  the  case  in  which  all  jobs 
can  be  completed  on  time.  From  the  due  date  constraints,  we  have  r}  +  pj  <  dj,  for  V;  6  J2- 
Hence, 

£  r,"  <  £  d3 ■,-  £  Pj  =  UB0 

J6J2  j€Ji  j€Ji 

However,  UB0  provides  a  weak  upper  bound.  Ahmadi  and  Bagchi  (1987)  derive  a  stronger 

bound  for  the  case  of  a  two-machine  flow  shop  in  which  all  job  due  dates  are  equal.    To 

our  knowledge,  this  is  the  only  other  upper  bound  available  currently  for  multiple  machine 

problems  addressing  similar  objectives.  We  will  denote  this  upper  bound  by  U B.  The  second 

set  of  experiments,  consequently,  addressed  a  2-machine  flow  shop  with  equal  job  due  dates 

and  compared  the  solution  value  yielded  by  MSP  with  respect  to  UB. 

5.1      Experimental  Design 

The  first  set  of  experiments  considered  a  5-machine  system  with  20  jobs.  Each  job  was 
assigned  5  operations,  and  the  machine  visitation  sequence  was  assigned  randomly  though 
successive  operations  of  a  given  job  were  processed  in  different  machines.  Operation  pro- 
cessing times  varied  uniformly  in  the  interval  (0,99).  The  due  date  of  dj  of  a  given  job  j  was 
determined  by 

dj  =  F(Y,Pj) 

j 

F  was  sampled  from  a  uniform  distribution  in  the  interval  (F  —  RF/2,  F  +  RF/2).  F  and  R 
respectively  control  the  tightness  and  the  variability  of  job  due  dates.  In  this  study,  two  sets 
of  due  date  tightness  levels  were  used.  Four  levels  of  tight  due  dates  were  obtained  by  setting 
F  at  0.05,  0.1,  0.15  and  0.20.  These  values  resulted  in  positive  tardiness  values  (although  at 
the  value  of  0.20,  the  tardiness  values  were  very  small).  Similarly,  loose  due  dates  at  four 
levels  were  obtained  at  F  values  of  0.30,  0.40,  0.50,  and  0.60.  At  each  of  these  eight  levels, 
two  values  of  R-O.b,  and  1.5,  were  used  to  provide  a  total  of  16  combinations  of  due  date  j 
tightness  and  variability. 

Five  instances  of  each  problem  scenario  were  randomly  generated.    Each  instance  was 
solved  using  two  approaches.  The  first  approach  employed  the  solution  procedure  discussed 
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in  Section  3  that  considered  both  primary  and  secondary  objectives.  The  second  approach 
considered  only  the  primary  objective  of  minimizing  total  tardiness.  The  solution  values  with 
respect  to  total  tardiness  and  the  sum  of  job  release  times  obtained  under  both  approaches 
were  recorded  and  averaged  over  the  five  problem  instances  for  reporting  purposes.  In  order 
to  restrict  the  computational  costs  within  reasonable  limits,  the  Modified  Operation  Due 
Date  (MOD)  rule  [see,  for  example,  Baker  (1984)]  was  used  for  solving  the  mean  tardiness 
problem  under  both  approaches.  In  total,  the  first  set  of  experiments  considered  80  problems. 
The  second  set  considered  a  2-machine  flow  shop.  The  smaller  problem  considered  25 
jobs  while  the  larger  problem  had  50  jobs.  The  operation  processing  times  were  selected 
from  a  uniform  distribution  in  the  interval  (0,99).  All  jobs  had  the  same  due  date  d  which 
was  determined  by 

j 
Five  values  of  F-0.6,  0.7,  0.8,  0.9,  and  1.0  were  used  to  generate  increasingly  loose  due 

dates.    [As  reported  in  Ahmadi  and  Bagchi's  (1987)  study  as  well,  we  found  that  F  <  0.6 

led  to  due  date  infeasibility  in  many  cases.]  As  in  the  first  set  of  experiments,  five  problem 

instances  were  solved  for  each  scenario.  For  each  instance,  the  ratio  of  the  sum  of  job  release 

times  obtained  by  solving  ICP  to  the  upper  bound  derived  by  using  Ahmadi  and  Bagchi's 

approach  was  recorded.    The  results  report  the  average  as  well  as  the  minimum  and  the 

maximum  values  of  these  ratios  over  the  five  problem  instances.    In  all,  the  second  set  of 

experiments  considered  50  problems. 

5.2      Experimental  Results 

Tables  1  and  2  give  the  results  of  the  first  experiment.  These  tables  compare  the  performance 
of  ICP  and  MOD  with  respect  to  total  tardiness  and  the  sum  of  job  release  times.  For  better 
comparison,  the  reported  total  tardiness  values  are  normalized  with  respect  to  the  sum  of 
job  processing  times.  Similarly,  the  sum  of  job  release  times  is  normalized  with  respect 
to  the  sum  of  job  due  dates.  In  both  tables,  z1?  and  z2  denote  normalized  total  tardiness 
and  normalized  sum  of  release  times  respectively.  Table  1  gives  the  results  for  the  case  of 
tight  due  dates,  while  Table  2  deals  with  the  loose  due  dates.  In  both  tables,  the  last  column 
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reports  the  average  solution  time  for  the  proposed  algorithm  on  an  IBM  3090-600  mainframe 
computer. 

Table  1  indicates  that  ICP  results  in  substantial  improvement  in  22  values  over  MOD  even 
when  due  dates  are  quite  tight.  This  improvement  generally  increases  as  due  dates  become 
progressively  looser.  Note  also  that  ICP  gives  marginally  better  total  tardiness  values  as 
well.  This  is  due  to  two  factors.  First,  because  of  the  job  indexing  scheme,  ICP  employs 
a  second  tie-breaking  rule.  Whenever  two  or  more  operations  are  found  to  have  the  same 
MOD  value  at  the  time  a  scheduling  decision  is  to  be  made,  it  selects  the  operation  for 
the  job  with  the  earlier  due  date.  Second,  note  that  ICP  first  schedules  jobs  in  J7i,  while 
MOD  considers  all  jobs  to  be  schedulable  at  any  given  time.  Therefore,  in  constructing  in 
nondelay  schedule,  MOD  frequently  takes  up  a  nonurgent  job  for  processing  (if  by  doing  so, 
machine  idleness  is  avoided).  This  may  delay  the  processing  of  an  urgent  job  that  arrives  at 
the  machine  soon  thereafter. 

It  can  be  seen  from  Table  2  that  relative  performance  of  ICP  over  MOD  is  superior  at 
higher  values  of  R.  This  is  due  to  the  fact  that,  with  greater  dispersion  in  job  due  dates, 
the  final  schedule  contains  a  larger  number  of  blocks.  Hence,  a  larger  portion  of  jobs  is 
completed  close  to  their  due  dates  with  consequent  increase  in  release  times  as  well. 

Table  3  presents  the  results  of  the  second  set  of  experiments.  The  values  reported  are 
the  ratios  of  the  ICP  solution  value  to  the  upper  bound  U B.  The  results  indicate  the 
effectiveness  of  the  algorithm  for  solving  ICP2,  at  least  in  the  case  of  a  2-machine  flow  shop. 
In  all  problem  instances,  the  solution  value  is  found  to  be  within  4.8%  of  the  upper  bound. 

6     SUMMARY  ; 

This  paper  examines  the  effectiveness  of  employing  input  control  as  a  mechanism  for  deferring 
job  release  in  a  just-in-time  system.  The  input  control  decision  is  modeled  as  a  dual  objective 
problem  of  minimizing  total  job  tardiness  and  maximizing  the  sum  of  job  release  times  in  a 
lexicographic  manner. 

The  proposed  solution  method  uses  a  constructive  approach  in  which  jobs  are  grouped 
into  blocks  for  simultaneous  processing.  The  membership  of  each  block  is  updated  in  an  iter- 
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ative  manner  by  solving  a  sequence  of  subproblems.  Our  computational  experience  indicates 
that  this  approach  yields  substantial  improvement  in  job  release  times  without  affecting  the 
tardiness  values  adversely.  This  result  is  of  interest  to  an  operating  manager  who  is  responsi- 
ble for  controlling  work-in-process  levels  and  earliness  costs  in  addition  to  meeting  due  dates 
effectively. 
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Appendix  1 

Notation 

j  Job  index,  j  =  1, . . . ,  N 

m  Machine  index,  m  =  1, . . . ,  M 

dj  Due  date  of  job  j,  j  =  1,. . .  ,N 

Pj  Processing  time  of  job  j,  j  =  1, . . . ,  N 

rj  Release  time  of  job  j,  j  =  1, . . . ,  N 

Sj  Start  time  of  job  j,  j '  =  1, . . .  ,iV 

Cj  Completion  time  of  job  j,  j  =  1, . . . ,  N 

Oj  Set  of  pairs  of  adjacent  operations  in  job  j,  j  =  1, . . . ,  N 

Nj  Number  of  operations  in  job  j,  j  =  1, . . . ,  N 

Tj  Tardiness  of  job  j  =  max(0,  Cj  —  dj),  j  =  1, . . .  ,iV 

Ej  Earliness  of  job  j  =  max(0,  d2  —  Cj),  j  =  1, . . . ,  N 

Pji  Processing  time  of  operation  i  in  job  j,  j  =  1, . . .  ,  N,  i  =  1, . . .  ,Nj 

rji  Release  time  of  operation  i  in  job  j,  j  =  1, . . . ,  N,  i  =  1, . . . ,  Nj 

Cji  Completion  time  of  operation  k  in  job  j,  j  =  1, . . . ,  N,  i  =  1, . . . ,  Nj 

Sji  Start  time  of  operation  k  in  job  j  j  =  1, . . . ,  N,  i  =  1, . . . ,  iVj 

fim  Set  of  all  operations  processed  on  machine  m,  m  =  1, . . . ,  M 
Ui  0  £  Vm   if  operation  i  in  job  j  requires  machine  m 
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Appendix  2 
Proofs 

Proof  of  Remark  1 

In  an  optimal  solution,  if  J2j€J  Tj  >  0,  then 

ii)  Tj  <  c(Ji\j)  +  Pj  -  dj,Vj  e  Ju   and 
Hi)  J\  starts  at  time  zero. 

PROOF:  i)  Clearly,  ^2j£j  >  YljeJi-   Hence,  if  the  condition  does  not  hold,  there  is  at  least 

one  j  E  Ji  such  that  Tj  >  0.  Let  u  be  the  job  with  the  earliest  start  time  in  J2-  Left  shifting 
all  j  and  all  jobs  that  start  before  j  such  that  u  starts  at  the  completion  time  of  J\  will 
not  increase  the  completion  time  of  these  jobs,  and  because  tardiness  is  a  regular  measure, 
will  not  increase  their  tardiness  as  well.  Clearly,  jobs  that  start  after  j  is  completed  will 
remain  unaffected.  Repeating  this  argument  for  all  jobs  with  positive  tardiness  provides  the 
required  result. 

ii)  Let  this  condition  not  hold  for  some  j  £  J .  Then  modifying  the  schedule  such  that 
j  starts  at  time  c(Jt\j)  will  not  increase  its  tardiness,  while  it  will  increase  its  release  time. 
Repeating  this  argument  as  often  as  required  gives  the  desired  result. 

Hi)  Follows  in  a  straightforward  manner  from  the  observation  that  the  tardiness  of  any 
job  cannot  increase  if  it  is  started  earlier. 

Proof  of  Remark  2 

In  each  block  Bk,  k  =  l,...,/£  in  an  optimal  schedule,  there  exists  one  critical  job  j  such 
that 

Cj  =  dj. 
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PROOF:    Because  each  job  j  in  Ji  is  early,  Cj  <  dj.    Suppose  the  condition  stated  in  the 
remark  is  not  satisfied  in  block  Bk-  Let  r3  =  d3  —  c7,  and 

u  =  arg  rnin]£Bk{Tj} 

Then  deferring  the  start  time  of  each  job  in  Bk  by  tu  increases  job  release  times  without 
delaying  any  job.  Such  a  right  shift  will  result  in  u  being  completed  exactly  at  its  due  date, 
and  it  will,  therefore,  be  the  critical  job. 
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Appendix  3 

Algorithm  for  Solving  BMRP 

1.  Initialization: 

a)  Renumber  all  jobs  as  per  EDD. 

b)  Set  Jx  =  0,  J2  =  J. 

c)  Determine  the  maximal  partition  of  J2  and  the  resulting  set  of  job  blocks  {  Bk  } ,  k  = 
1,...,A'. 

2.  Sequencing  Jobs  within  Blocks: 

a)  For  each  block  Bk  £  J2,  solve  BMRP.  Set 

Dk  =  min3€Bk{d3\c3  -  Ck)i    ?k  =  ™axjeBk{c3}  -  ™>injeBk{sj] 

If  Dk  >  Pk,  go  to  Step  3.  Otherwise,  go  to  Step  2b. 

b)  Set  J\  =  J\  U  Bk-  Solve  ICP1  for  Jx-  Update  Pk  based  on  the  revised  schedule. 
If  Dk  <  Pk,  then  set  J2  =  J2\Bk. 

3.  Block  Sequencing:  If  Jx  ^  0,  set  Rx  =  0.  Starting  with  the  first  block  in  J2,  schedule 
each  block  Bk  €  J2  such  that  Ck  =  Dk-  If  Sk  =  Ck  —  Pk  <  Ck-x  for  any  block  k,  go  to 
Step  4.  Otherwise,  stop. 

4.  Block  Merging:  Set  K  =  K  —  1.  If  k  =  1,  merge  Bx  with  Jx,  reschedule  jobs  within  Jx 
by  solving  ICP1.  Else,  merge  Bk  with  Bk-x-  Solve  BMRP  for  the  block  thus  formed. 
Update  Pk-x  and  Dk-x  accordingly.  Go  to  Step  3. 

Algorithm  for  Obtaining  the  Maximal  Partition 

In  the  following  algorithm,  it  is  assumed  that  all  jobs  are  numbered  in  the  EDD  order.  Also, 
J  denotes  the  set  of  jobs  already  assigned  to  blocks,  and  Jc  =  J  \  J.  The  algorithm  steps 
are  given  below. 
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1.  Set  k  =  1;  J  =  0;  Jc  =  J\  j\  =  1.  Compute  d3  —  pj,  j  =  1, . . .  ,  iV. 

2.  Find  the  largest  job  index  j2  in  ^c  such  that 

<*j"2  ~  Ph  <  dh  • 
Assign  jobs  {ji, . . . ,  j2]  to  block  k. 

3.  Set  Jc  =  Jc  \  {ju. . .  ,j2).    If  Jc  =  0,  set  A'  =  k  and  stop.    Else,  set  k  *—  k  +  1, 
ii  < —  J 2  +  1  •  Go  to  Step  2. 
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TABLE  1 

Comparison  of  ICP  and  MOD  -  Values  of  (z\,  z2) 
Tight  Due  Dates 


R 

F 

(z\\z2) 

z2{ICP)/z2(MOD) 

CPU 

MOD 

ICP 

0.5 

0.05 

1.76;  0.54 

1.72;  0.71 

1.31 

2.442 

0.10 

0.63;  0.34 

0.62;  0.42 

1.24 

1.785 

0.15 

0.12;  0.26 

0.12;  0.34 

1.31 

2.002 

0.20 

0.00;  0.22 

0.00;  0.35 

1.59 

3.804 

1.5 

0.05 

1.52;  0.54 

1.49;  0.70 

1.30 

1.757 

0.10 

0.51;  0.39 

0.50;  0.46 

1.18 

1.851 

0.15 

0.08;  0.31 

0.08;  0.39 

1.26 

2.217 

0.20 

0.02;  0.24 

0.02;  0.62 

2.58 

4.952 

TABLE  2 

Comparison  of  ICP  and  MOD  -  Values  of  (z\,  z2) 
Loose  Due  Dates 


R 

F 

izi'i  zi) 

z2(ICP)/z2{MOD) 

CPU 

MOD 

ICP 

0.5 
1.5 

0.30 
0.40 
0.50 
0.60 
0.30 
0.40 
0.50 
0.60 

0.00;  0.16 
0.00;  0.12 
0.00;  0.10 
0.00;  0.08 
0.00;  0.12 
0.00;  0.12 
0.00;  0.10 
0.00;  0.08 

0.00;  0.95 
0.00;  0.77 
0.00;  0.83 
0.00;  0.88 
0.00;  0.74 
0.00;  0.82 
0.00;  0.88 
0.00;  0.91 

5.94 

6.42 

8.30 

11.00 

6.17 

6.83 

8.80 

11.34 

1.889 
1.831 
1.654 
1.204 
1.359 
0.912 
1.758 
0.066 
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TABLE  3 


Values  of  Z2{ICP)/UB  for  a  Two  Machine  Flow  Shop 


NJ 

F 

Minimum 

Average 

Maximum 

25 

0.6 

0.952 

0.973 

0.996 

0.7 

0.963 

0.979 

0.997 

0.8 

0.970 

0.983 

0.997 

0.9 

0.975 

0.986 

0.998 

1.0 

0.978 

0.988 

0.998 

50 

0.6 

0.956 

0.972 

0.992 

0.7 

0.965 

0.978 

0.994 

0.8 

0.972 

0.982 

0.995 

0.9 

0.976 

0.985 

0.996 

1.0 

0.979 

0.987 

0.997 
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Figure  1  -  Equivalence  of  BMRP  and  CTP 
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